Thermally-induced magnetic phases in an Ising spin Kondo lattice model 

on a kagome lattice at 1/3-filling 
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Numerical investigation on the thermodynamic properties of an Ising spin Kondo lattice model on 
a kagome lattice is reported. By using Monte Carlo simulation, we investigated the magnetic phases 
at 1/3-filling. We identified two successive transitions from high-temperature paramagnetic state to 
a Kosterlitz-Thouless-like phase in an intermediate temperature range and to a partially disordered 
phase at a lower temperature. The partially disordered state is characterized by coexistence of 
antiferromagnetic hexagons and paramagnetic sites with period \/3 x v3. We compare the results 
with those for the triangular lattice case. 
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I. INTRODUCTION 



Geometrical frustration has been one of the major top- 
ics in the condensed matter physics for a long time. A 
variaty of peculiar phenomena arc emergent from the 
macroscopic degeneracy of the ground state and its lift- 
ing by perturbative interactions |l| . Thermal or quantum 
fluctuation also plays a significant role on the system, and 
in some cases, it gives rise to an exotic phase or unusual 
order. A typical example is found in the Ising antifcrro- 
magncts on a triangular lattice. When the interaction is 
restricted to the nearest neighbor spins, the geometrical 
frustration prevents the system from long-range ordering, 
and the system remains disordered down to zero temper- 
ature with extensive degeneracy and associated residual 
entropy 0-[H . The extensive degeneracy can be lifted by, 
e.g., further neighbor interactions and external magnetic 
field, leading to a variety of magnetic phases. 

An example of such unusual phases is the partially dis- 
ordered (PD) state. The PD state has a peculiar mag- 
netic order in which paramagnetic sites coexist with mag- 
netically ordered sites. A possibility of such state was 
first pointed out in the triangular-lattice Ising model in 
the presence of second-neighbor ferromagnetic interac- 
tion in addition to the nearest-neighbor antiferromag- 
netic interaction. Within a mean-field calculation, a 
three-sublattice PD state was obtained between high- 
temperature paramagnetic state and low-temperature 
three-sublattice ferrimagnetic state Q . In this PD state, 
the antiferromagnetically ordered moments on the honey- 
comb subnetwork coexist with the paramagnetic ones on 
the remaining sites. However, subsequent Monte Carlo 
(MC) calculations have indicated that the PD state is 
unstable in the purely two-dimensional systems, and is 
taken over by a Kosterlitz-Thouless (KT) phase 0^]. 
The situation was essentially the same with additional 
further-neighbor interactions @ and also on a kagome 



(a) 



o o 




'Electronic address: ishizuka@aion.t.u-tokyo. ac.jp 



FIG. 1: (Color online) Schematic pictures of the magnetic 
structure of PD states with (a) q = and (b) y/3 x y/3 or- 
dering. The shaded region in (a) shows the unit cell for the 
kagome lattice, while the shade in (b) is the magnetic super- 
cell for the \fZ x y/3 order. The numbers in (a) denote v in 
eq. ([2]), and A,B,C in (b) denote a in eq. ©. 



lattice [l0( . These results imply that the PD state is hard 
to be stabilized in two-dimensional Ising spin models. 

On the other hand, recently, the authors have reported 
stabilization of a PD state in two dimensions by intro- 
ducing a coupling to itinerant electrons [Tl| . By a nu- 
merical investigation using MC simulation technique for 
an Ising spin Kondo lattice model on a triangular lat- 
tice, they showed that a three-sublattice PD state ap- 
pears with opening a (pseudo) charge gap at finite tem- 
peratures above a phase-separated region. The results 
indicated that the non-perturbativc effect of spin-charge 
coupling plays crucial roles in stabilizing PD. 

To address the generality of such PD state stabilized by 
spin-charge coupling, we here conduct numerical investi- 
gation of an Ising spin Kondo lattice model on a kagome 
lattice. In the localized spin systems, kagome Ising anti- 
ferromagnets have several distinct features compared to 
the triangular ones. When the exchange interaction is 
restricted to the nearest-neighbor antiferromagnetic one, 
the kagome model also exhibits macroscopic degeneracy 
in the ground state. However, the residual entropy in 



1 



the kagome case is S = 0.502 [l4[, which is considerably use the pseudo-spin Q defined for each unit cell 
larger than the value in the triangular case, S = 0.323 Q. 
This is due to the stronger frustration in the kagome lat- 
tice consisting of corner-sharing triangles, compared to 
the edge-sharing triangular lattice; indeed, the spin cor- 
relation in the ground state shows an exponential decay 
in the former, while it is a power-law decay in the lat- 
15| . For this reason, it is interesting to make a com- 
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parative study between the kagome and triangular cases 
in the presence of the spin-charge coupling. It will pro- 
vide a good test for examining the stability of PD state. 

In this contribution, we investigate the thermodynamic 
behavior of an Ising spin Kondo lattice model on a 
kagome lattice by using a MC simulation technique. As 
the result, we confirmed two successive transitions: one 
is from paramagnetic phase to a KT-like magnetic phase, 
and the other is from the KT-like state to v3 x a/3 PD 
state [Fig. |Tfb)] . The results together with the previous 
ones for the triangular lattice model suggest that the sta- 
bilization of PD states by spin-charge coupling is rather 
general in two-dimensional Ising spin Kondo lattice mod- 
els. 



II. MODEL AND METHOD 

We consider a single-band Kondo lattice model on a 
kagome lattice (see Fig. [1} with localized Ising spin mo- 
ments. The Hamiltonian is given by 



H 



H.c. 
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The first term represents hopping of itinerant electrons, 
where Ci a (c\ a ) is the annihilation (creation) operator of 
an itinerant electron with spin a =^, 4- at ith site, and t is 
the transfer integral. The sum is taken over nearest 
neighbor sites on the kagome lattice. The second term 
is the onsite interaction between localized spins and itin- 
erant electrons, where of represents the ^-component of 
itinerant electron spin and Si = ±1 denotes the localized 
Ising spin at ith site; J is the coupling constant (the sign 
of J does not matter in the present model). Hereafter, 
we take t = 1 as the unit of energy, the lattice constant 
o = l, and the Boltzmann constant ks = 1- 

To investigate thermodynamic properties of this 
model, we adopted a MC simulation which is widely 
used for similar models [l6| . The calculations were con- 
ducted up to the system sizes with the number of the 
unit cell iV = 6x6, 6x9, and 9x9, under the pe- 
riodic boundary conditions (the total number of lattice 
sites is given by N s = 3N). Thermal averages of phys- 
ical quantities were calculated for typically 4300-9800 
MC steps after 1700-5000 MC steps of thermalization. 
In the following, we focus on the electron density at 
» = J2i,A4*Ci*)/N 8 = 1/3. 

To distinguish the PD and other magnetic phases, we 



where S]^ is the Ising spin on ^th site in mth unit cell 
[Fig. (Ha)]. Then, we define the sublatticc pseudo-spin 
moment by 
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where a = A, B, C denotes three sublattices for the 
pseudo-spins [Fig. QJb)]. The PD states shown in Fig. [1] 
are characterized by a finite M° = (M^M^Mf) par- 
allel to (s/3/2, I/n/2,0), (0, a/2,0), or their threefold 
symmetric directions around the z-axis. On the other 
hand, the ferrimagnetic states, in which the paramag- 
netic spins in the PD states are aligned fcrromagnctically, 
are distinguished by a finite M° along ( a/2/3, v2, 1 / \/3) , 

(2^/2/3,0, — I/a/3), or their threefold symmetric direc- 
tions around the z-axis. Hence, to parametrize them, we 
use the azimuth parameter defined by 



i> a = (M^) 3 cos6<^ 



(4) 



where 0m° is the azimuth of M a in the xy-plane, and 
M x « = 3(Af^) 2 /8 [M« = {(Al?) 2 + (Af«) 2 } 1/2 ] [ll],[l7|. 
The parameter ip a takes a negative value and ip a — ¥ 
—27/64 for the ideal PD ordering, while it becomes pos- 
itive and tp a — > 1 for the ferrimagnetic ordering; ip a = 
for both paramagnetic and KT phases in the thermody- 
namic limit N — > oo. 

For both q = and s/3 x s/3 orders shown in Fig. [TJ 
M" y , \M"\, and V" are independent of a. In the follow- 
ing calculations, we found that all the quantities are in- 
dependent of a; hence, we show the averaged values over 
a > °' = E Q O q /3 (O a = M* y , |Af z Q |, V Q )- To distinguish 
the q = order, we also calculate the net pseudo-spin 
moment 



M 



(5) 



and corresponding M xy , \M Z \, and ip. In addition, we 
compute the susceptibilities for the pseudo-spin moments 
for each case to detect the phase transitions. 



III. RESULTS AND DISCUSSION 

Figurc[2]shows the results of MC calculation at n = 1/3 
and J = 2. As shown in Fig.[2][a), M' monotonically in- 
creases as decreasing temperature, and shows a rapid in- 
crease around Tkt = 0.056(3). In addition, it exhibits a 
small shoulder around T C PD = 0.014(2) before approach- 
ing the value at the lowest temperature. The two anoma- 
lies are more clearly observed in the corresponding sus- 
ceptibility x'xy plotted in Fig. &b) ; x' xy shows a divergent 
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FIG. 2: (Color online) MC results for (a) M' xy and \M' Z \, (b) 
Xxy an d x'z! an d (c) ip' . The data are for J — 2 and with the 
system sizes JV = 6 x 6, 6 x 9, and 9 x 9 at n = 1/3. 



peak at Tkt and a hump structure at T PD . The results 
imply two successive transitions at Tkt and T PD . 

First, we consider the temperature region T < T C PD . 
At the lowest temperature of our calculation, M' ap- 
proaches y/2, while \M' Z \ is essentially zero in the ther- 
modynamic limit, as shown in Fig. [2ja). In addition, 
ijj' shows a sharp decrease in this region, from ip' = to 
~ —0.4, as shown in Fig.[2Jc). The nonzero ip' indicates a 
spontaneous breaking of six-fold rotational symmetry of 
M"; specifically, the negative value approaching —27/64 
suggests that the system exhibits PD. Correspondingly, 
as shown in Fig. [U[b), x'z shows a monotonic increase as 
decreasing temperature, which is ascribed to the pres- 



ence of paramagnetic spins in the PD state. In contrast 
to M 1 , no increase in the net moment M xy nor M z is 
seen in the entire range of calculation; ip is also zero. We 
also found that the spin structure factor exhibits a peak 
corresponding to the v3x v3 order (not shown). Hence, 
we conclude that the system exhibits the PD state with 
period a/3 x y/3 for T < T PD . 

Next, we examine the intermediate temperature region 
T PD < T < Tkt- As shown in Fig. Efb), x' xy shows a 
divergent peak at Tkt corresponding to a rapid rise of 
M' xy in Fig. HJa), which is a clear indication of a phase 
transition. M' xyl however, exhibits a considerable finite 
size effect in this temperature region r PD < T < Tkt- 
On the other hand, \M' Z \ and xp' shows almost no change. 
In particular, ip' is extrapolated to zero within statistical 
errors in the limit of N — ¥ oo, indicating that the inter- 
mediate region shows no rotational symmetry breaking 
in respect to M°. Similar behaviors were observed in 
the KT phase with quasi long-range order in the Isinj 
antiferromagnets on triangular and kagome lattices 
IToj and in the Kondo lattice model on a triangular lat 
tice ll| . Hence, we consider that the intermediate phase 
for T PD < T < Tkt is of KT type. 

Consequently, our MC data indicate that the system 
exhibits two successive phase transitions in the calculated 
temperature range: one is from high-temperature param- 
agnetic phase to the intermediate KT-like phase, and the 
other is from the KT-like phase to the low-temperature 
PD state with period V3 x -\/3- One point to be ad- 
dressed here is an anticipated phase transition from the 
PD state to the true ground state. The PD state re- 
tains residual entropy of ~ ^ log 2 associated with the 
1/3 paramagnetic moments. Hence, it is unlikely for the 
PD state to be the ground state of the present model be- 
cause the degeneracy will be lifted by long-range RKKY 
interactions induced by the spin-charge coupling [l8l - [20j . 
However, in our MC simulation at n = 1/3, there is no 
indication of further phase transition from PD down to 
T ~ 10 -3 t. This implies that the energy scale of the 
relevant RKKY interaction is extremely small and that 
the true ground state is nearly degenerate with the PD 
state. Similar feature was also reported recently in an 
Ising spin Kondo lattice model on a triangular lattice, 
where the PD state survives down to an extremely low 
temperature at n = 1/3 fill ]. 

Finally, let us compare the present results with those 
for the model on a triangular lattice [ll[ ■ In the kagome 
lattice case, the transition temperature to the PD state, 
T C PD = 0.014(2), is roughly ten times smaller than that 
in the triangular lattice case, T PD = 0.130(4). This is 
presumably because of the stronger frustration in the 
kagome lattice case. Namely, the antiferromagnetic or- 
der develops on the honeycomb network in the case of 
PD on the triangular lattice case, whereas it appears on 
the disconnected hexagons in the kagome lattice case, as 
shown in Fig.[2Jb); the former can be accommodated by 
nearest neighbor interactions, but the latter needs fur- 
ther neighbor interactions. Hence, the difference of T PD 
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can be understood since the RKKY interaction becomes 
weaker for further neighbors in general. In the kagome 
lattice case, instead, the KT-like phase appears in a wide 
temperature range above the PD phase, whereas a direct 
transition from paramagnetic to PD phase was observed 
for the triangular lattice model at 1/3-hlling. This might 
also be ascribed to the difference of frustration; higher 
entropy due to the stronger frustration may contribute 
to the prevalence of the KT-like phase. 

IV. CONCLUSIONS 

To summarize, we have investigated the magnetic 
properties of an Ising spin Kondo lattice model on a 
kagome lattice at 1/3-filling. To elucidate the thermo- 
dynamic properties of this model, we employed an unbi- 
ased Monte Carlo simulation and calculated the pseudo- 
spin moments, corresponding susceptibilities, and az- 
imuth parameter. All the parameters calculated consis- 



tently suggested that there are two successive phase tran- 
sitions, from the high-temperature paramagnetic phase 
to an intermediate-temperature Kostcrlitz-Thouless-like 
phase, and to the low-temperature -s/3 x y/3 partially- 
disordered state. Comparison between the present re- 
sults and those for the triangular lattice case implies that 
the partial disorder is stabilized in a wide class of two- 
dimensional Kondo lattice models under frustration. 
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